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We analyse matrix elements of AB = 2 operators by combining QCD results with the ones obtained in the 
static limit of HQET. The matching of all the QCD operators to HQET is made at NLO order. To do that we 
have to include the anomalous dimension matrix up to two loops, both in QCD and HQET, and the one loop 
matching for all the AB = 2 operators. The matrix elements of these operators are relevant for the prediction of 
the B — B mixing, Bs meson width difference and supersymmetric effects in AB = 2 transitions. 



1. INTRODUCTION 

Matrix elements of AB = 2 operators cru- 
cially enter theoretical determinations of impor- 
tant phenomenological quantities. Due to the 
high mass of the b quark and to the present com- 
puting power, it is impossible to simulate b quarks 
directly on the lattice. Two approaches have been 
employed to overcome this problem: 

• Perform relativistic simulations with sev- 
eral masses around the c-quark mass and 
extrapolate the results to physical 6-quark 
mass. 

• By means of an effective theory, the heavy 
degrees of freedoms can be integrated out 
and only the light modes are simulated. 
So far, HQET and NRQCD have been 
employed. In HQET an expansion in 
^QCD is made whereas in NRQCD the 
limit Tj ^ 1 is taking, being v the heavy 
quark velocity. 

In this work we have combined, for the first 
time, full QCD and HQET simulations to com- 
pute the matrix elements of 4-quark AB — 2 op- 
erators. To do so appropriately, the computation 
of the matching coefficients of the operators be- 
tween QCD and HQET at next-to-leading order 
(NLO) accuracy is needed. 

*Talk given by Juan Reyes. 



We have considered the following basis of 4- 
quark AB — 2 operators: 
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with i, j colour indices, and q stands for either d- 
or s- light quark flavour. The first of the above 
operators enters the Standard Model (SM) de- 
scription of the B^ — B° mixing amplitude [Q, 
whereas O2 and O3 are relevant for the SM esti- 
mates of the relative width difference in the neu- 
tral S-meson system, (Ar/r)^ The matrix 
elements of all operators parametrize supersym- 
metric effects in AB = 2 transitions [^. 

It is usual to express the matrix elements of 
the operators (|l]) in terms of the so-called B- 
parameters, which are introduced as a measure 
of the deviation from the vacuum saturation ap- 
proximation (VSA), namely 

(B0|Oi(A.)|i?,") - b,mljlB,{p), 

= hxmlJlB,{^l), (2) 
2 < i < 5 

with b = {8/3, -5/3, 1/3, 2, 2/3} and x = 
m?g / {mb{fi) + mq{fi))^ . The hat symbol denotes 
operators renormalizcd in some renormalization 
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scheme at the scale /z. To determine the values 
of the B-parameters i?i_5(/z), we have combined 
the results of a (quenched) QCD numerical simu- 
lation on the lattice with 3 values of the heavy 
quark mass, in the range of heavy-light pseu- 
doscalar masses mp G (1.7, 2.4) GeV, with HQET 
results in order to constrain the extrapolation to- 
wards the physical point, ms^/^- This extrapola- 
tion is guided by the Heavy Quark Scaling Laws 
(HQSL) which are the propertly ones of HQET. 
Therefore, to use the HQSL we have to match the 
S-parameters obtained in QCD onto the HQET 
ones. We refer the reader to ref. [|j where this 
issue is explained in great detail. 

2. SIMULATION DETAILS 

We enumerate, briefly, the main elements of our 
lattice simulations: 

• For full QCD: the details can be found in 
refs. |^-||. The simulation is performed in 
a lattice of the size 24^ x 48, at /3 = 6.2 
corresponding to = 2.7 GeV, with 
the non-perturbatively improved Wilson ac- 
tion. The number of gauge configurations 
is 200, with 3 values of the heavy and 
3 values of the light quark masses, cor- 
responding to the Wilson hopping param- 
eters: Kq e {0.1344,0,1349,0.1352}, and 
KQ G {0.125,0,122,0.119}. 

• For HQET: the details are in ref. The 
simulation is performed in a lattice of the 
size 24'^ X 40, aX j3 — 6.0 corresponding 
to a^^ — 2.0 GeV, with the tree level 
Clover improved Wilson action. The num- 
ber of gauge configurations is 600, with 
3 values of the light quark masses, corre- 
sponding to the Wilson hopping parameters 
Kq e {0.1425,0.1432,0.1440}. 



3. EXTRAPOLATION 
B-MESON 



TO 



THE 



Once the B-parameters of the renormalized op- 
erators, in QCD and HQET, are obtained from 
lattice simulations in some scheme and at some 
scale, the matching between the two theories 



reads: 



W, 



C{mp) ■ WHQET[mp,fA'^ ■ B{fi) 
1 \ 



O 



(3) 



where W^qc£)[m27 /^i]^^ is the matrix encoding 
the full QCD evolution from a scale /zi to ^2 
of all five AB = 2 B-parameters , likewise for 
^HQBT 1/^2, Ml] in HQET. The matrix C is the 
matching matrix of the B-parameters between 
QCD and HQET. These matrices are specified in 
^ . B{mp,^) is a five-component vector which 
collects all five AB = 2 B-parameters renor- 
malized at the scale /i simulated with a heavy 
quark mass corresponding to a pseudescalar me- 
son mass of mp. -B(/i) is the corresponding vector 
in HQET. 

In order to account for the logarithm depen- 
dence in eq. (H) we put the HQET evolution 
matrix in the l.h.s. and construct the quantity 
$(top,^): 

<l>(mp,/i) = 
■B{mp,y) 

which can be fit either freely as 



(4) 



^{mp,fi) = ao(/i) 



mp 



(5) 



where ao(/i) and di{fi) are the fit parameters, or 
by constraining it by the static HQET results, i.e. 



<i>{mp,n) = a^(^) 



mp 



(6) 



where the coefficient cl'qIp) is completely con- 
strained by the static value, B{^), so that one 
can probe the term 0{l/m^). In figure 1 we 
show the two extrapolation of eqs. |^ and ^ for 
the first three B-parameters. In the analogous 
plots for <I>4 and <I>5 can be found. As a result, 
we obtain the HQET values of the _B-parameters, 
i.e. ^(mp^^^, /i), which are then to be matched 
back onto their QCD counterparts. The final re- 
sults in the constrained case are presented in table 



3 



1.2 
1.0 
0.8 
0.6 

1.2 
1.0 
0.8 
0.6 

1.2 
1.0 
0.8 
0.6 



■ ■ ■ - -^^ --^ f-# 



O <I>^ (nip, m^) 



I I I I 



I I I I 



I I I I 



_1 — — — ■ — 



~$ ^ — ^ 



• BJm.) 



I I I I 



I I I I 



I I I I 



r 

_L_ 



O (nip, m^) 



0.0 



0.5 



1.0 



1.5 



1/M^ 



Figure 1. Extrapolation to the physical Bd meson 
mass (squared symbols) in the inverse heavy meson 
mass. The dotted line corresponds to the uncon- 
strained linear extrapolation for each of the compo- 
nents "l>i(mp,m(,) from our data (empty circles) to 
^i(rnB^,m^,) (empty square). The result of the con- 
strained extrapolation (filled squares) by the static 
HQET B-parameters (filled circles) is marked by the 
dashed line. 



1, where the first error is the statistical one and 
the second the systematic error, due to the un- 
certainty in the lattice renormaUzation constants, 
both in QCD and HQET (see ref. % for more de- 
tails). 

This work has been supported by the European 
Community's Human potential programme under 
HPRN-CT-2000-00145 Hadrons/LatticeQCD. 



Table 1 

B-parameters defined in eq. renormalized at fi = 
rub — 4.6 GeV and in two renormaUzation schemes: 
RI/MOM and MS of ref. 

Scheme 





RI/MOM 


MS(NDR) [|10| 


B'f\m,) 


0.87(4) (tl) 


0.87(4) (tl) 


B^im,) 


0.82(3)(4) 


0.79(2)(4) 


Bf(m,) 


l.G2(6)(9) 


0.92(6)(8) 


B^im,) 


1.16(3) (t?) 


1.15(3) (t'r) 


Bf\m,) 


1.91(4) (If) 


1.72(4) (tf ) 


B['\m,) 


0.86(2) (tl) 


0.87(2) (tl) 


Bi'\m,) 


0.83(2)(4) 


0.80(1)(4) 


Bi'\m,) 


1.03(4)(9) 


0.93(3)(8) 


Bi'\m,) 


1.17(2) (t?) 


1.16(2) (t?) 


Bi'\m,) 


1.94(3) (if) 


1.75(3) (if) 
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